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Kinetic equations for magnetic nano particles dispersed in a viscous liquid are developed and 
analyzed numerically. Depending on the amplitude of an applied oscillatory magnetic field the par¬ 
ticles orient their time averaged anisotropy axis perpendicular to the applied held for low magnetic 
held amplitudes and nearly parallel to the direction of the held for high amplitudes. The transition 
between these regions takes place in a narrow held interval. In the low held region the magnetic 
moment is locked to some crystal axis and the energy absorption in an oscillatory driving held is 
dominated by viscous losses associated with particle rotation in the liquid. In the opposite limit the 
magnetic moment rotates within the particle while its easy axis being nearly parallel to the external 
held direction oscillates. The kinetic equations are generalized to include thermal huctuations. This 
leads to a signihcant increase of the power absorption in the low and intermediate held region with 
a pronounced absorption peak as function of particle size. In the high held region, on the other 
hand, the inclusion of thermal huctuations reduces the power absorption. The illustrative numerical 
calculations presented are performed for magnetic parameters typical for iron oxide. 

PACS numbers: 75.40.Gb, 75.40.Mg, 75.75.Jn, 75.60. Jk 


I. INTRODUCTION 

In recent years magnetic nano particles are discussed 
intensively in connection with biomedical applications, 
in particular magnetic hyperthermia. The use of these 
particles in hyperthermia is based on the fact that power 
is absorbed in an applied oscillating magnetic field. The 
goal is to achieve a high absorption rate with particle 
parameters suitable for biomedical applications. This 
restricts the amplitude of the applied magnetic field to 
100 — 200 Oe and its frequency to 50 — 1000 kHz [Il-Q. 

Theoretically the dynamics of magnetic nano particles 
has been studied in recent years in very many papers for 
the case that the particles are fixed in space. The dy¬ 
namics is then reduced to that of magnetic moments in 
external fields for which the stochastic Landau-Lifshitz- 
Gilbert (LLG) equation introduced by Brown Q is of¬ 
ten used. This approach has gained increased interest 
recently because of its application to single-domain par¬ 
ticles. In these nano particles (NPs) the internal degrees 
of freedom can be neglected at low temperatures making 
it possible to describe the whole particle as one macro 
spin of constant length Q. Its dynamics is expected to 
be well described by a classical approach. The dynam¬ 
ics of these macro spins is greatly influenced by thermal 
fluctuations originated from coupling of the spins to the 
surrounding medium (the heat bath). The theory of ther¬ 
mal fluctuations of small magnetic particles is one of the 
fundamental issues of modern magnetism [^. 

As far as magnetic hyperthermia is concerned a re¬ 
cent comprehensive study is available which focusses on 
the problem of obtaining NPs with optimal parameters 
with respect to high energy absorption rates needed in 
biomedical applications Q- However, the theoretical ap¬ 
proach used in this work relies on various simplifications 
like linear response theory, two level approximation for 


magnetic NPs with uniaxial anisotropy or the concept of 
a frequency dependent coercivity. A more fundamental 
study based on the LLG equations has also been pub¬ 
lished recently Q which avoids the approximations men¬ 
tioned. In both of these works the NPs are assumed to 
be fixed in space. 

In biomedical applications, however, a new aspect 
should be considered, the possibility for magnetic par¬ 
ticles to move in the viscous liquid they are immersed 
in. The behavior of such particles has been studied in 
some detail, in particular in the theory of magnetic flu¬ 
ids iOl. In these systems the magnetic moment may 
rotate within the particle with respect to some crystal 
axis -i.e. Neel relaxation - and rotate together with the 
particle with respect to the liquid. Both processes con¬ 
tribute to the power absorption of NPs when placed in a 
time varying magnetic field. Theoretical studies have as¬ 
sumed for the most part [1, that both processes occur 
separately leading to the concept of effective relaxation 
times treated within linear response theory [HHil- This 
heuristic approach was critically discussed in detail in [l^ 
because it certainly oversimplifies the physical situation, 
since it does not take into account the rich dynamics of 
magnetic NPs. Moreover the field amplitudes interesting 
for applications are beyond the linear response regime 
making a detailed investigation of the dynamics of mo¬ 
bile NPs in finite alternating fields necessary. 

Therefore, in the present paper this dynamics is stud¬ 
ied starting from the well established LLG equations for 
magnetic NPs suitable modified for taking the rotational 
degrees of freedom into account. Our aim is to establish 
a set of kinetic equations for studying the dynamics of 
mobile single domain NPs at finite temperatures. A sim¬ 
ilar approach has been discussed and applied to magnetic 
hyperthermia very recently in 0111. In these papers, 
however, sets of kinetic equations are proposed which 
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both violate in limiting cases basic conservation laws to 
be valid on physical grounds making the results obtained 
doubtful. An important achievement of the present inves¬ 
tigation is the proposal of a set of kinetic equations which 
are in accordance with necessary conservation laws. 

The NPs to be considered are assumed to have an uni¬ 
axial anisotropy characterized by an unit vector n par¬ 
allel to the easy axis. Their magnetic moment has a 
fixed absolute value and is oriented parallel to the 
unit vector e. Numerical solutions of the kinetic equa¬ 
tions proposed show that depending on the amplitude of 
an applied oscillating magnetic field basically two differ¬ 
ent regions occur in the stationary state. In a low field 
region n and e are nearly parallel and their time aver¬ 
age is perpendicular to the direction of the applied field. 
In the high field region, on the other hand, n is nearly 
parallel to the direction of the applied field while e os¬ 
cillated between the states ±z where z is the direction 
of the applied field. Such a scenario has been found and 
discussed already in M- The results obtained, however, 
differ from our findings in many details. In particular, 
we found an extremely sharp transition between these 
regions at an applied field amplitude of about one half 
of the coercive field accompanied by a dramatic increase 
of the area of the hysteresis loop when passing from the 
low field to the high field region. 

Thermal fluctuations play an important role in the dy¬ 
namics of magnetic NPs. The kinetic equations proposed 
are generalized accordingly. Numerical solutions of these 
equations are focused on a calculation of the power ab¬ 
sorbed by the system, particularly on its dependence on 
the size of the particles. For low field amplitudes a reso¬ 
nance type absorption peak as function of particle size is 
found at finite temperatures. The corresponding power 
absorption is much higher than without thermal fluctu¬ 
ations. This enlarged power absorption originates from 
temperature driven switching processes of the magneti¬ 
zation relative to the NPs and therefore is strongly de¬ 
pendent on the volume of the NPs. For increasing field 
amplitudes these absorption peaks first broaden and then 
go over to a monotonous behavior in the high field region. 
In this region thermal fluctuation reduce the power ab¬ 
sorption especially for very small particle sizes. Results 
are compared with those obtained for an ensemble of NPs 
with fixed randomly distributed anisotropy axis. In this 
case the power absorption is generally smaller than in 
the case of mobile particles. In view of applications in 
hyperthermia the illustrative numerical calculations are 
performed for magnetic parameters typical for iron oxides 
at room temperature. 


II. BASIC EQUATIONS 


A nano particle considered in the present paper is mod¬ 
eled as a uniformly magnetized spherical rigid body (RB) 
with uniaxial magnetic anisotropy so that its Hamilto¬ 


nian is given by 

■H = l/2eu}^ - D{n-nf. ( 1 ) 


0 denotes the moment of inertia of the RB, u its angu¬ 
lar velocity, the magnetic moment with magnitude fig, 
an anisotropy energy, B the external magnetic field 
and n an unit vector parallel to the anisotropy axis of the 
RB. This vector often called director is firmly attached 
to the RB so that its equation of motion is given by 

9n 

-=u,x„. (2) 

The differentiation in Eq. m is performed in a coordi¬ 
nate system fixed in space, the laboratory frame. 

The equations of motion for the angular momentum 
L = 0 u; of the RB and for the spin momentum S asso¬ 
ciated with the magnetic moment are determined by 
the corresponding torques given by 


dt 


= M, 


(3) 


and 



(4) 


The torque M; can be obtained from the energy change 
following an infinitesimal rotation of the anisotropy axis 
n, i5n = X n, according to M/ • S(f = —SH with the 
result that [l^ 


M; = 2 • n) n X /X. (5) 

The torque Mg can be obtained in a similar way from 
the energy change following an infinitesimal rotation of 
the moment, 5fj, = 6(px fi, with the result that 

Mg = fi X B + 2 D{fi ■ n) fi X n. (6) 


The sum of both is the net torque in which the intrin¬ 
sic contributions stemming from the anisotropy energy 
cancel each other. The equation of motion for the total 
angular momentum S -|- L is thus given by 


5(S-bL) 

dt 


fjb xB 


(7) 


showing that it is conserved in zero field as it should be. 
Note that angular momentum L and spin momentum S 
are coupled through the anisotropy energy so 

that the conservation of the total angular momentum is 
nontrivial. 

So far we have only considered the conservative part 
of the equations of motion. Dissipative contributions 
are added phenomenologically. Here we have to distin¬ 
guish between internal and external relaxation processes. 
External processes arise due to the interaction of the 
NPs with the surrounding fluid in the form of rotational 
damping and thermal fluctuations. 
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Internal relaxation processes to be considered first oc¬ 
cur due to the motion of the spin relative to the NP. For 
this spin dynamics we introduce a Gilbert damping term 
proportional to S x Because it arises from the mo¬ 
tion of the spin S it must vanish if this vector does not 
move relative to the RB. Therefore in the present case the 
time derivative in this term must be interpreted as time 
derivative in the moving frame fixed to the RB. Thus in 
terms of derivatives in the laboratory frame the Gilbert 
damping term should be written as -d (S x (^ — a; x S)) 
where a denotes a damping constant. This term has to 
be added to the right hand side of Eq. (g]) which together 
with Eq. (6) can be rewritten as 

— =/I X Be-5(S X (—-a; X S)) (8) 

where Be denotes an effective field, 

Be = B -I- 2 D{iJ, ■ n) n. (9) 


The damping term discussed is a torque acting on the 
spin moment due to the interaction of the magnetic mo¬ 
ment with the other degrees of freedom of the RB. There¬ 
fore, the opposite momentum is transferred to the rota¬ 
tional degrees of freedom, i.e. the Gilbert damping term 
has to be subtracted from the right hand side of Eq. (5) 
resulting in 

— = n) fj, X n +a{S X {— - u} X S)). (10) 

It follows that Eq. (7) is still valid so that for an iso¬ 
lated particle the total angular momentum L -|- S again 
is conserved in zero field as it should be irrespectively of 
internal damping processes. The results obtained are in 
contrast to Ref. [ij in which the second term on the 
right hand side of Eq. (10) is missing. They also dif¬ 
fer from equations proposed i n IlSlI which were already 
criticized of being incorrect in [14| . 

Further sources of damping are external relaxation pro¬ 
cesses stemming from the rotational motion of the NP in 
the surrounding liquid. We assume a damping term lin¬ 
ear in u} 


N; = (11) 

with friction coefficient ^ = GrjVd- Here, ry is the dynami¬ 
cal viscosity while Vd denotes the hydrodynamic or total 
volume of the particle, i.e. it is assumed that the parti¬ 
cle consists of a magnetic core region of volume Vm with 
radius eventually covered by a nonmagnetic surfac¬ 
tant layer. Here it is assumed that the surfactant layer 
is firmly attached to the magnetic core. 

The torque Ni has to be added to the right hand site 
of Eq. (10) so that the equation of motion for L now 
reads 

f)T _ og 

— = — 2 D{fi-n) /XX n+a(S X xS)) — ^u}. (12) 


The equation of motion for the total angular momentum 
S + L is therefore given by 


9(S + L) 
dt 


= /X X B -I- N/. 


(13) 


To proceed we first solve Eq. (8) for ^ with the result 

-<5(S X (/X X Be)-t-S X w)). (14) 

An explicit equation of motion for L can then be obtained 
by inserting Eq. (14) into Eq. (13). 

Finally we introduce an unit vector e in the direction of 
the magnetic moment, e = /x//is .The spin momentum 
S can be expressed as /x = —7 S where 7 denotes the 
gyromagnetic ratio. We express S and /x in terms of e, 
define D = Hs^D, a = a VS • S and end up with the 
following set of equations 

— = — -{e X Be + ae X (e X Be) 

ot 1 -I- V 

— ex (u; X e) -I - e x u: (15) 

7 79 


dt ~ j dt 


HsG X B — ^o). 


(16) 


These equations together with Eq. (I2|) are the basis of the 
following calculations. At low temperatures they specify 
the dynamics of the NP completely. At elevated tem¬ 
peratures thermal fluctuations have to be added to these 
equations. This will be discussed in section IV. 

Note that the quantity a appearing in Eq. (15) is the 
dimensionless damping parameter usually used in the lit¬ 
erature [l6l - [l^ . In the present paper the kinetic equa¬ 
tions proposed will be discussed exclusively in the low 
frequency limit because only this limit is relevant for 
biomedical applications. Note, however, that the equa¬ 
tions proposed are not restricted to this limit. 


III. APPLICATION TO A NP IN A VISCOUS 
LIQUID 

A cartesian coordinate system spanned by unit vectors 
X, y and z and fixed in space is used. The applied alter¬ 
nating magnetic field is assumed to be aligned parallel to 
one of these axis, the z-axis, and is given by 

B = Ho cos(w t) z (17) 

where w = 27r/ is the angular frequency of the exter¬ 
nal driving field. For practical applications in hyper¬ 
thermia the frequency / is of the order of 50 — 1000 
kHz. For the numerical calculations it is convenient to 
use the dimensionless field amplitude bg = BojBc where 
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Be = 2D/Us denotes the coercive field. We scale the 
time t according to t = tuji where w; = 2'yD/^s is the 
Larmor frequency and we introduce the reduced angular 
velocity ft = w/w;. The resulting dimensionless equa¬ 
tions are easily seen to depend only on the parameters 
R = /^isBc, S = 6r]^Vdyn/^J's and F = 27r//wi. 

To perform illustrative numerical calculations we 
choose magnetic parameters typical of the particles of 
iron oxides (magnetite) [T3 |. The anisotropy energy D 
can be expressed as D = Ki ■ Vm and Ki is the magnetic 
anisotropy constant given by Ki = 10“^ J/m^ for mag¬ 
netite. The magnetic moment fig can be expressed as 
/is = Mg ■ Vm with Mg = 4 • lt)^A/m. With these values 
one obtains for the coercive field = 50 mT and for 
the reduced friction coefficient S = 2.64- 10^ • •q/'q^uater ■ 
Vd/Vm with q^ater = 10“^- kg/ms. 

In the following we use 

q — q/qwater * VdjVm ( 1 ^) 

because the dynamical viscosity q and the hydrodynamic 
volume Vd enter in this combination throughout. Note 
that q is the sole parameter depending on the interaction 
of the NP with the surrounding fluid in the absence of 
thermal fluctuations. 

For the parameters underlying our numerical calcula¬ 
tions the kinetic equations can be simplified if one takes 
into account [ij, [l^ that because of the presence of a 
viscous damping together with the smallness of the NP 
it is possible to neglect the term involving the moment 
of inertia, assuming R = 0, so that u) can explicitly be 
obtained from Eq. (HU). To study the accuracy of this ap¬ 
proach we performed calculations with and without using 
this approximation and found no detectable differences in 
the results. 

The kinetic equations are solved numerically. Initial 
conditions for e, n and at time t = 0 have to be spec¬ 
ified. Usually the reduced magnetization is chosen to be 
e = X while the easy axis of the NP is set parallel to x-l-z 
and 0 = 0. The system of equations is damped so that 
a stationary state is reached finally. 

Under the influence of an oscillating driving field the 
magnetization switches initially on very small time scales 
of the order of Tinit ^ 10 — 20 towards the direction of 
the easy axis into a state parallel to the effective held 
Bg. Note that on these short time scales the applied 
held is practically constant. After reaching this position 
magnetization and anisotropy axis start oscillating with 
the frequency of the applied field accompanied by a slow 
relaxation into a stationary state. This relaxation process 
is shown for two values of the amplitude of the driving 
field in Fig. 1. Here, Uzm denotes averaged values of the 
z-component of the anisotropy axis averaged over periods 
of the applied field for / = 50 kHz. The symbols in 
Fig. 1 represent these quantities n^m which are set at 
odd multiples of to/ 2 where tq denotes the period of the 
driving field which in the present case is given by tq = 
1.76 TO®,!, e. huge as compared to . Closed symbols 
refer to bg = 0.54 while open symbols refer to the smaller 



FIG. 1. (color online) Averaged z-component of the 
anisotropy axis for bg = 0.54 (upper three curves) and 
bo = 0.48 (lower three curves). ?) = 10 (circles, black), 30 
(squares, red) and 60 (diamonds, blue). / = 50 kHz, a = 0.5. 

field bg = 0.48. Results for three different values of q are 
shown: circles refer to q = 10, squares refer to ?) = 30 
while diamonds refer to q = 60 . 

Fig. 1 shows that the averaged component rizm decays 
to zero for bg = 0.48 which means that in the average n 
orients perpendicular to the driving field in the stationary 
state. On the other hand for bg = 0.54 in the average 
n orients nearly parallel to the direction of the driving 
field. An increase in the viscosity only slows down this 
relaxation into the stationary states. 

The results discussed are typical in the sense that in¬ 
dependently of details in the parameters two states are 
found, a state for small fields in which the direction of 
the anisotropy axis averaged over field periods is orien¬ 
tated perpendicular to the applied field and a state for 
larger fields in which it is oriented nearly parallel to the 
field. We now discuss more details of the dynamics in 
the stationary states. 

In Fig. HI components of e and n are shown in the 
stationary state as function of reduced time r for field 
amplitude bg = 0.52 in the upper panel and for bg = 0.53 
in the lower panel. Shown are the amplitudes after about 
20 field cycles sufficient for reaching stationarity (for 
q = 10). For bg = 0.52 the x-components of both the 
anisotropy axis, rij, (dashed green), and the reduced mag¬ 
netization, Cx (solid black), are close to one with small 
oscillations in time. As was discussed above in this case 
the NP is orientated in such a way that its easy axis is on 
the average parallel to the x-axis, i.e. perpendicular to 
the driving field. The magnetization is practically locked 
to the NP. This can be seen even more clearly by consider¬ 
ing the scalar product e-n (not shown) between the scaled 
magnetization and the director of the NP which in the 
present case turns out to deviate from one by less than 5 
%. The components (dotted dashed blue) and Cz (dot- 
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FIG. 2. (color online) Cx (solid black), Cz (dotted red), Ux 
(dashed green), (dotted-dashed blue), a = 0.5, fj = 10, 

/ = 50 kHz. Upper panel: bo — 0.52, lower panel bo = 0.53. 


ted red), however, both oscillate around zero while the 
corresponding y-components are extremely small. Thus, 
with numerical accuracy n lies in the x — z plane with 
a large component parallel to the x-axis and a smaller 
^-component which oscillates around zero. Therefore the 
angular velocity O (not shown) is practically parallel to 
the ?/—direction. In the present case its y-component 
turns out to be of the order of 10“^ while the other com¬ 
ponents of H are negligible. 

For larger fields a completely different behavior is ob¬ 
served. The lower panel in Fig. [5]shows that for fog = 0.53 
the z-component of n (dashed dotted blue) is nearly 
aligned with the direction of the driving field (the z-axis) 
which means that the particle has rotated into a direc¬ 
tion with the easy axis nearly parallel to the driving field. 
The magnetization e, on the other hand, rotates between 
the states ±z in the x — z-plane. It is important to note 
that in both cases the amplitude of the driving field is 
well below the coercive field be = 1.0. 



FIG. 3. (color online) Field dependence of averaged com¬ 
ponents: Uzm (dots, black), Cxm (squares, red), (n1)m (dia¬ 
monds, blue), a = 0.5, / = 50 kHz, rj = 10. 

The two cases discussed are typical in the sense that 
there exists a low field region in which e is locked within 
the RB with an easy axis oriented in the average perpen¬ 
dicular to the driving field and a high field region in which 
e rotates within the RB. For a quantitative analysis we 
calculated the time averaged x-component of the reduced 
magnetization in the stationary state, as well as the 
time averaged z-component of n averaged over periods of 
the oscillating field as a function of the amplitude of the 
driving field. From the discussion above we expect e^m 
close to one for small fields and zero in the large field re¬ 
gion while the opposite behavior is expected for Uz- The 
results shown in Fig. [3] confirm this. A remarkable sharp 
transition is observed between these regions at a thresh¬ 
old field bth around 0.525. Also shown in Fig. |3]is (nl)m 
which jumps from a small value below the threshold to a 
large value above it showing that in the low-field region 
there are finite oscillations of Uz around its mean value 
nzm = 0.0 while in the high field region these oscillations 
are around a finite mean value. 

The switching behavior of immobile NPs is fre quen tly 
discussed within the Stoner - Wohlfarth model in 
which the energetically stable and metastable positions 
of the magnetization of a single domain NP with uniax¬ 
ial anisotropy in the presence of a static magnetic field 
are considered. It turns out that the minimum static 
field needed to reverse the magnetization at zero temper¬ 
ature depends on the direction of the applied field with 
respect to the easy axis as described by the well-known 
Stoner-Wohlfarth astroid. The smallest switching field is 
given by i?c/2 and therefore close to the threshold field 
found in the present case. Note, however, that a com¬ 
plete agreement with this result is not expected since on 
the one hand time-dependent magnetic fields as consid¬ 
ered here reduce the required minimum field below the 
Stoner-Wohlfarth limit [2ll - l^ and on the other hand the 
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FIG. 4. (color online) Hysteresis loops for / = 50 kHz 
(upper panel) and for / = 200 kHz (lower panel). With 
increasing order of the loop area field valnes are bo = 
0.3,0,52,0,53,0.56,0.8. 


mobile NPs are not forced into a specific direction with 
respect to the applied field. Nevertheless the resemblance 
of minimum switching field and threshold field seems to 
be remarkable. 

Additional calculations show that on the one hand 
side the clear separation between a low field region with 
nzm = 0 and a high field region with a finite value of Uzm 
is observed for all values of the parameters studied. On 
the other hand the precise value of the threshold field sep¬ 
arating these regions although always found to be around 
bo ^ 0.5 turned out to depend on parameters like damp¬ 
ing, viscosity or frequency. Thus the driven dynamical 
system can settle in various stationary states which be¬ 
come accessible by changing the system parameters. A 
detailed study of this complex dynamics is not followed 
any further here because we are mainly interested in the 
finite temperature behavior of the system for which these 
details are expected to be not so important. 


The transition between the two regions discussed above 
is accompanied by a dramatic change in form and area 
of the hysteresis loops. Typical loops are shown in Fig. 
4 for low frequencies, / = 50 kHz (upper panel) and 
/ = 200 kHz (lower panel). The loop for bo — 0.3 (black) 
has an ellipsoidal shape with a small area. It widens 
slightly when going up with the field to 6o = 0-52 (green). 
Increasing the field further to 6 = 0.53 a drastic increase 
of the loop area is observed together with a change to a 
more rectangular shape. With increasing frequencies the 
loop area decreases significantly in the low field region 
while in the high field region no dramatic changes are 
found as can be seen in Fig. S) 

The area of the hysteresis loops is proportional to the 
power absorbed. Our calculations show that a dramatic 
increase of the absorption within a narrow field interval 
occurs when passing from the low to the high field region 
around bo ^ 0.5. In the high held region the magne¬ 
tization switches relative to the NP in contrast to the 
low held region where only oscillations of the NP rela¬ 
tive to the viscous huid occur. Thus the large increase of 
the energy absorption observed originates obviously from 
the transition of oscillatory dynamics to switching of the 
magnetization relative to the NP. We note in passing that 
due to the reduced squareness of the hysteresis loop an 
estimate of its area only from data of the dynamical coer- 
civity as was done in Ref. Q is not meaningful for small 
held amplitudes. 

The absorbed power can be obtained either by measur¬ 
ing the loop area or alternatively from 0, [Tj, 

QTD 

P = -M«(e—). (19) 

The use of this expression is more convenient especially 
at hnite temperature where large thermal huctuations 
have to be taken care of. The brackets in Eq. (IT^ de¬ 
note a time average over full held cycles of the driving 
held in the stationary state. In the present case the re¬ 
duced averaged absorbed power pav = V/2KiVm can be 
expressed explicitly as 

Pav = tbbo {ezSin{Cbt)). (20) 

Note that the quantity e and therefore Pav are indepen¬ 
dent of the volume of the NPs in the limit i? —>■ 0 as a 
closer inspection of Eqs. (15) and (16) shows. 

In Fig. [5]we show Pav with the same sets of parameters 
as used above for different values of the driving frequency 
/. The solid symbols show results for a reduced viscosity 
t) = 1.0 while open symbols represent results for ?) = 50. 
The expected dramatic increase of the absorbed power 
when passing through the threshold region is clearly vis¬ 
ible. Furthermore a strong increase with increasing fre¬ 
quency is observed. On the other hand, an increase of 
the effective viscosity fj leads to a significant reduction 
of Pav especially for high frequencies and in the low field 
region where Pav is dramatically reduced down to about 
6 • 10^ for ij = 50 for all frequencies shown. The rea¬ 
son for this finding obviously is the reduced mobility of 
the NPs leading to less frictional losses. 
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FIG. 5. (color online) Averaged absorbed power pav vs. bo 
for different values of /. Closed symbols are for p — 1.0, open 
symbols for fj = 50. / = 1000 kHz: diamonds, green; / = 500 
kHz: squares, red; / = 100 kHz: circles, black, a = 0.2. 

As was noted above for a fixed value of the effective vis¬ 
cosity fj the reduced absorbed power pav is independent of 
the magnetic volume of the NPs. Thermal fluctuations, 
on the other hand, are strongly volume dependent lead¬ 
ing therefore not only to quantitative but also to quali¬ 
tative changes of the results presented. In particular, we 
will show in the next sections that a resonance type ab¬ 
sorption peak as function of particle size occurs for small 
amplitudes of the driving field. It originates from tem¬ 
perature driven switching processes of the magnetization 
which have a strong influence on the power absorption. 

IV. THERMAL FLUCTUATIONS 

Thermal fluctuations are expected to have a signifi¬ 
cant impact on the behavior of NPs, especially if their 
radius becomes very small. Here we assume that the spin 
dynamics is governed by the stochastic Landau-Lifshitz- 
Gilbert (LLG) equation in which a fluctuating field C 
is added to the effective field Bg. This fluctuating field is 
generally assumed to be Gaussian distributed with zero 
mean and correlator 

( 0 ( 0 )Cm(t))c = Si,mS{t)2akBT/{p,-f). (21) 

In the same spirit the equation of motion for a; has to 
be modified by adding a fluctuating torque e to the right 
hand site of Eq. (16) which leads to Brownian rotation 
of the RB in the liquid [l3, [111 . e is also assumed to be 
Gaussian distributed with zero mean and with correlator 

m 

{eii0)em{t))c = Si,mS{t)2^kBT. ( 22 ) 

In both of these equations I and m label cartesian com¬ 
ponents of the fluctuating fields. 


In the following the stochastic equations are solved nu¬ 
merically with emphasis on a calculation of the power 
absorption by a set of NPs. Methods well known from 
literature are used [IMl- The power absorption is cal¬ 
culated using Eq. (fTOl) in which the brackets now denote 
an average over a very large number of field cycles in 
order to reduce the effect of thermal fluctuations. Note 
that alternatively one can split this into an ensemble av¬ 
erage over a set of NPs with a correspondingly reduced 
number of field cycles. In the following we average over 
about 4000 field cycles altogether after reaching a sta¬ 
tionary state. We found that this value is sufficient in 
general to obtain reproducible results. 

Thermal fluctuations depend strongly on the volume 
of the NP. It is therefore of interest to study the - 
dependence of Pav in detail. A reduction to dimensionless 
units reveals that for a fixed set of magnetic parameters 
Ki and Mg the kinetic equations depend only on two 
parameters, the reduced viscosity fj introduced in Eq. (1181) 
and rewritten as 

fj — V/ Vwater ‘ {j’d/'^mj • 

and a scaled temperature factor, 

T = Tlrl. (24) 

In these equations denotes the radius of the magnetic 
core of the NP and its hydrodynamic radius. 

These relations show in particular that an increase of 
the magnetic core region reduces the temperature fluctu¬ 
ation while an increase of the hydrodynamic volume of 
the NP has a similar effect as an increase of the viscos¬ 
ity (for fixed r^). It also means, for instance, that for 
a given ratio rd/vm the parameter T goes to zero for a 
large particle size Vm so that Pav becomes independent 
of Tm and approaches the T = 0 value. Eor intermediate 
particle sizes, on the other hand, a complex dependence 
on Tm and 6o is expected. 

A. Small fields 

To elucidate these dependences we first show in Pig. [S] 
data for pav obtained in the low field region for 6o = 0.1 
and for 6o = 0.4 (inset of Fig. [6]). The horizontal bro¬ 
ken lines show Pav for T — 0 which is independent of 
magnetic radius as was discussed before. For finite 
temperatures, however, sharp resonance type absorption 
peaks are observed exceeding the T = 0-values signifi¬ 
cantly. The peaks for T = 100 K (circles, red) move to 
larger values of when increasing the temperature to 
T = 300 K (squares, blue). For increasing the absorp¬ 
tion approaches the T = 0-values as is expected from the 
general discussion. 

The reason for these observations is a temperature 
driven switching of the magnetization leading to an in¬ 
crease of the absorption in a restricted r^-region. For 
very small switching is very fast resulting in a nearly 
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FIG. 6. (color online) Averaged absorbed power vs. magnetic 
radius Vm for bo = 0.1. The inset shows results for bo = 0.4. 
For both graphs: T = 0 K (horizontal lines, black), T = 100 
K ( circles, red), and T = 300 K (squares, blue); f] = 1.0, 
a = 0.2, / = 500 kHz. 


reversible hysteresis curve with negligible energy absorp¬ 
tion. For larger values of Neel -type relaxation sets 
in leading to large absorption rates. For even higher 
values of this type of relaxation is blocked result¬ 
ing in a strong reduction of the absorption which, for 
very large values of r^, approaches the T = 0 - value 
as expected. The shift of the resonance-type absorption 
curve to higher values of with increasing temperature 
confirms this picture. 

When increasing the field amplitude from bo = 0.1 
to bo = 0.4 (inset of Fig. 6) the resonance type peaks 
broaden considerably and show not only an increase of 
the maximal absorption by a factor of more than 10 but 
also a shift of its position to higher values of with 
increasing value of the field amplitude. 

For a deeper insight into these mechanism we calcu¬ 
lated for bo = 0.1 the scalar product e • n in the sta¬ 
tionary state. While for such small driving fields this 
quantity is very close to one at T = 0 independently 
of Tm it oscillates at finite temperatures irregularly be¬ 
tween the states ^ ±1 on very small time scales. These 
temperature driven switching processes are very fast for 
small values of and slow down for increasing parti¬ 
cle size up to a point where switching is blocked. This 
slowing down is shown in Fig. 7 for two values of 
around the maximum of the absorption curve. The data 
are obtained in the stationary state over a time inter¬ 
val of length At ^ 16 tq where tq denotes the period 
of the driving field. The scalar product e • n is aver¬ 
aged over very small time intervals of length S = tq/IOO, 
s = (e • n)^, just in order to smoothen the data some¬ 
what for drawing purposes. This quantity s is plotted 
in Fig. 7 together with the driving field which is shown 
as wavy lines. It is seen that switching of the magne¬ 


FIG. 7. (color online) Averaged scalar product s (see text 
for definition) for two values of vs. r over time intervals 
of length Ar together with driving field (wavy line, blue), 
bo = 0.1, j) = 1.0, a = 0.2, / = 500 kHz. 

tization relative to the easy axis of the NP slows down 
significantly when increasing the size of the NP. 

The energy absorption is correlated with these temper¬ 
ature driven switching processes. For very small values 
of Tjn these processes are very fast and the response of 
the magnetization with respect to the driving field is al¬ 
most reversible leading to small power absorption. Large 
power absorption is observed when the frequency of the 
external fields roughly matches the switching frequency 
which takes place in a narrow - interval. A further 
increase of the particle size leads to a blocking of the tem¬ 
perature assisted switching resulting in a sharp decrease 
of the power absorption towards its T = 0 - value. 

1. Immobilized NPs 

It is interesting to compare the results obtained with 
those for an ensemble of immobilized magnetic NPs. Out 
of the large number of theoretical papers on the proper¬ 
ties of such systems some of them have been devoted 
^ecifically to the problem of magnetic hyperthermia 
H i Hi Hi- It is therefore of relevance to compare 
our results with those obtained for an ensemble of im¬ 
mobilized NPs. To this purpose we present numerical 
calculations on the power absorption based on our ki¬ 
netic equations. They are reduced to the case of immo¬ 
bilized NPs using the same set of material parameters 
as before. The ensemble consists of 150 NPs with fixed 
anisotropy axis randomly distributed in space. In Fig. 8 
Pav is shown for such an ensemble of NPs for bo = 0.1 
and otherwise the same set of parameters as used in Fig. 
6 . 

A comparison of the results shows that the mobile NPs 
have a significant higher absorption rate throughout. On 
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FIG. 8. (color online) Averaged absorbed power for an ensem¬ 
ble of immobile NPs with randomly distributed anisotropy 
cixis (circles, black) and for mobile NPs (squares, blue, same 
as in Fig. 6) for bo = 0.1. Inset shows results for b — 0.4. 
rj = 1.0, a = 0.2, / = 500 kHz. 


FIG. 9. (color online) Averaged absorbed power vs. magnetic 
radius Vm for different values of frequency. / = 1000 kHz ( 
circles, black), / = 500 kHz ( squares, red), / = 100 kHZ 
(diamonds, blue) for bo ~ 0.1, T = 300 K, rj — 1.0, a = 0.2. 
Horizontal broken line (black) shows pav for T = 0 K. The 
inset shows results obtained from Eq. 28. 


the other hand the similarity in these curves suggests 
that in both cases the main absorption mechanism are 
similar. This confirms our previous findings that a large 
power absorption is certainly due to magnetic switching 
processes of the magnetic moment relative to the NP. 
These switching processes are increasingly blocked for 
increasing r^. This leads to a dramatic decrease of the 
absorption for immobile NPs with increasing Vm- For 
mobile NPs, on the other hand, freezing of the switching 
processes with increasing also takes place but in this 
case Pav levels out at the constant T = 0-value the ori¬ 
gin of which is a significant absorption due to Brownian 
rotational motion. 


2. Phenomenological approach 


In [Hill a phenomenological approach was presented 
for heating magnetic fluids with alternating magnetic 
fields. It uses a simple relaxation equation for the mag¬ 
netization of a NP with an effective relaxation rate Ts 
consisting of two parts, Debye relaxation with relaxation 
time [13 


tb = 


ksT 


and Neel relaxation with relaxation time 


tn = \/kBT/D tt/loi eyi-p{D/kBT). 


(25) 


(26) 


Arguing that Brownian and Neel processes take place 
in parallel the effective relaxation time Ts in the linear 
response regime is supposed to be given by Q 

( 27 ) 


leading to a power dissipation of the following form 


Pav 


1 -I- (QTsY ' 


(28) 


An expression of this form has widely been used, see for 
instance EMI. This concept of an effective relaxation 
time, however, was criticized to be not strictly justified 
El- We note in passing that the energy absorption de¬ 
rived for a set of rigid NPs in linear approximation m 
also leads to an expression similar to that given in Eq. 


For a comparison of these heuristic results with results 
obtained from our kinetic equations we show in Fig. 9 
Pav for different frequencies in the limit of small field 
amplitudes, bo = 0.1. The inset of Fig. 9 shows pav as 
given by Eq. (28) for the same set of parameters demon¬ 
strating that this simple approximation leads indeed to a 
resonance-type absorption peak at a value of rm which, 
however, is only about half of what we found in our cal¬ 
culations. Moreover, the absorption goes to zero for large 
r^n contrary to our results where it approaches the con¬ 
stant T = 0 -values as was discussed above. Addition¬ 
ally, frequency dependence and width of the absorption 
peaks exceed those found in our approach. It is impor¬ 
tant to note that for this comparison no parameters are 
needed other than those we used in our numerical calcu¬ 
lations. The qualitative agreement of these results with 
those from our more fundamental approach is remark¬ 
able, in particular because there are no adjustable param¬ 
eters involved except for the over-all amplitude of Pav ■ It 
would be interesting to see under what assumptions and 
approximations our rather complicated equations can be 
reduced to such a simplified approach. 
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FIG. 10. (color online) Averaged absorbed power vs. mag¬ 
netic radius Tm for mobile NPs (squares, green) and for an 
ensemble of immobile NPs (circles, black). Dashed horizon¬ 
tal lines denote the corresponding T = 0 -values, feo = 0.65, 
T = 300 K, a = 0.2, r) = 2.0, / = 1000 kHz. 

B. Intermediate to high fields 

The broadening of the absorption peak and the shift of 
the position of its maximum to higher values of when 
increasing the field amplitude (see Fig. 6) continues for 
fields up to about 6o ~ 0.5. For larger fields, on the other 
hand, the peak structure disappears and a monotonic in¬ 
crease of Pav as function of Tm is observed for both mobile 
and immobile particles. This is demonstrated in Fig. 10 
for a field amplitude of 6o = 0.65 and driving frequency 
/ = 1000 kHz. The absorption increases as function of 
Tm and approaches for large Vm the T — 0-values from 
below. The power absorption pav for an ensemble of mo¬ 
bile particles is larger by more than a factor of two as 
compared to immobile particles. 

To be more specific, switching of the magnetization 
depends on temperature as expressed for instance in the 
Neel relaxation formula, Eq. (26), and on the value of 
the applied field, see for instance ref. [2^. As was dis¬ 
cussed above, for small fields Neel relaxation is important 
leading to an extremely strong dependence of the relax¬ 
ation time on the magnetic volume of the particle and 
on temperature leading in a restricted Tm - interval even 
to an increase of the power absorption above its T = 0 - 
value. For larger fields above 6o ~ 0.5, on the other hand, 
switching occurs already without thermal fluctuations at 
T — Q, c.f. Fig. 2. In this case thermal fluctuations re¬ 
duce the absorbed power for all values of in contrast 
to low fields where they are essential for high power ab¬ 
sorption. The largest reduction of Pav occurs for small 
values of where thermal fluctuations are most impor¬ 
tant. 

More detailed results for pav in the region of small 
are presented in Fig. 11, both for mobile and immobile 



r 

m 



FIG. 11. (color online) Upper panel: Averaged absorbed 
power vs. Vm for different field amplitudes. 

Lower panel: Averaged absorbed power vs. for an ensem¬ 
ble of NPs with fixed randomly distributed anisotropy axis. 
For both cases: bo = 0.65: triangles down, maroon; bo — 0.53: 
triangles up, green; bo = 0.4: squares, blue; bo = 0.3: dia¬ 
monds, red; bo — 0.2: circles, black. T = 300 K, a = 0.2, 
/ = 1000 kHz. 


particles. The ensemble of immobilized magnetic NPs 
consists of about 150 NPs with fixed anisotropy axis ran¬ 
domly distributed in space. All the other parameters are 
the same in both parts of Fig. 11. 

The results show that the mobile NPs have a signif¬ 
icant higher absorption rate throughout. On the other 
hand the similarity of these curves again suggests that in 
both cases the main absorption mechanism are similar. 
This confirms our previous findings that a large power 
absorption is certainly due to magnetic switching pro¬ 
cesses of the magnetic moment relative to the NP while 
viscous losses are relatively small except in those situa¬ 
tions where magnetic switching processes are blocked so 
that viscous losses are the only sources for power absorp- 
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tion. 

The most interesting result being relevant for appli¬ 
cations in hyperthermia, however, is the development 
of a very broad maximum of Pav as function of parti¬ 
cle size already for moderate field amplitudes of around 
bo ^ 0.4. The maximal value of Pav for bo = 0.4 is 
reached for ^ 12 nm and exceeds the T = 0-value 
by many orders of magnitude. Thus this temperature 
driven resonance type behavior has the important prac¬ 
tical consequence that large absorption rates are already 
obtained for relatively small fields and for values of 
out of a broad interval. These findings are of particular 
interest because in applications a distribution in particle 
sizes is expected so that absorption rates which are not 
too sensitive to particle sizes are advantageous. On the 
other hand, for fields above the threshold region around 
bo ~ 0.5 a monotonous increase of Pav as function of 
is observed so that NPs with large magnetic radius Vm 
are advantageous for power absorption. 


V. CONCLUSION 

In the present paper we proposed a set of equations 
for the dynamics of magnetic nano particles dissolved in 
a liquid. These equations are constructed on the basis 
of the LLG equation for classical spins together with a 
classical equation of motion for a particle in a viscous 
medium taking the rotational degrees of freedom into 
account. The obtained kinetic equations are such that 
physically necessary conservation laws are fulfilled in con¬ 
trast to other published work. 

Numerical calculations reveal clearly that without 
thermal fluctuation two distinct regions can be identified 
in a stationary state: for low amplitudes of the driving 
field the easy axis of the NP orients perpendicular to 
the field with a magnetization nearly fixed with respect 
to the easy axis of the NP. For larger fields the NP ori¬ 
ents in such a way that its easy axis is nearly parallel to 
the direction of the driving field while its magnetization 
switches relative to the NP along its easy axis. A very 
sharp transition between these regions is observed at zero 
temperature. 

The energy absorption calculated at finite tempera¬ 
tures increases considerably in the low field region and 
especially when going from the low field to the high field 
region. Sharp resonance type absorption peaks as func¬ 
tion of Tm are observed for small field amplitudes which 
broaden and finally turn over to monotonic behavior 
when increasing the field amplitudes. These results con¬ 
tradict the conclusions drawn in [l^ that in contrast to 
immobilized NPs the power absorption of NPs in a liquid 
increases monotonically with increasing particle diame¬ 
ter. Such a behavior is only observed for large enough 
field amplitudes. For smaller fields the choice for an op¬ 
timal particle size leading to a large power absorption 
is more delicate because large absorption rates are only 
obtained in rather narrow - intervals. The illustra¬ 


tive numerical calculations at finite temperatures were 
performed for parameters relevant for hyperthermia and 
can help finding optimal parameters for an effective en¬ 
ergy absorption by magnetic NPs. 

The present work is motivated by applications in hy¬ 
perthermia where magnetic NPs are likely to be free to 
rotate in a liquid. However, particles may also be im¬ 
mobilized at walls or inside biological cells [l^. Results 
we obtained for ensembles of NPs with fixed anisotropy 
axis show that as far as power absorption is concerned 
the difference between these two cases is not as impor¬ 
tant as expected before [l^, at least not for parameters 
interesting for hyperthermia. 

In this contribution we only considered single nano par¬ 
ticles embedded in a viscous liquid. For an ensemble of 
interacting NPs a generalization of the present formalism 
is obvious although the numerical effort increases consid¬ 
erably. In this case not only the interaction of the par¬ 
ticles has to be taken into account but also translational 
degrees of freedom because the particles will attract or re¬ 
pel each other affecting these degrees of freedom. Impor¬ 
tant in this context is the dipolar interaction which may 
lead to a clustering of particles depending on tempera¬ 
ture [1^. Additionally, in an alternating field switching 
processes of the magnetization may result in changes of 
the distance between particles leading to oscillations in 
the particle positions which may generate high frequency 
oscillations in the liquid [2^ . These interesting questions 
arising from particle interactions are left for future re¬ 
search. 
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